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Abstract 

The gravitational waveform (GWF) generated by inspiralling compact binaries moving in quasi- 
circular orbits is computed at the third post-Newtonian (3PN) approximation to general relativity. 
Our motivation is two-fold: (i) To provide accurate templates for the data analysis of gravitational 
wave inspiral signals in laser interferometric detectors; (ii) To provide the associated spin-weighted 
spherical harmonic decomposition to facilitate comparison and match of the high post-Newtonian 
prediction for the inspiral waveform to the numerically-generated waveforms for the merger and 
ringdown. This extension of the GWF by half a PN order (with respect to previous work at 2.5PN 
order) is based on the algorithm of the multipolar post-Minkowskian formalism, and mandates the 
computation of the relations between the radiative, canonical and source multipole moments for 
general sources at 3PN order. We also obtain the 3PN extension of the source multipole moments 
in the case of compact binaries, and compute the contributions of hereditary terms (tails, tails-of- 
tails and memory integrals) up to 3PN order. The end results are given for both the complete plus 
and cross polarizations and the separate spin-weighted spherical harmonic modes. 
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I. INTRODUCTION 



One of the most important sources of gravitational radiation for the laser interferometric 



detectors LIGO, VIRGO 



, l2j and LISA |3j is the inspiralling and merging compact binary 
system. Until the late inspiral, prior to merger, the gravitational waves are accurately de- 
scribed by the post-Newtonian (PN) approximation to general relativity [4J, while the late 
inspiral and subsequent merger and ringdown phases are computed by a full-fledged numer- 
ical integration of the Einstein field equations . A new field has emerged recently 
consisting of high-accuracy comparisons between the PN predictions and the numerically- 
generated waveforms. Such comparisons and matching to the PN results have proved cur- 
rently to be very successful jj], 



10. 



11 



121 ] . They clearly show the need to include high 



PN corrections not only for the evolution of the binary's orbital phase but also for the 
modulation of the gravitational amplitude. 

The aim of this paper is to compute the full gravitational waveform generated by inspi- 
ralling compact binaries moving in quasi-circular orbits at the third post-Newtonian (3PN) 
order 1 . By the full waveform (FWF) at a certain PN order, we mean the waveform in- 
cluding all higher-order amplitude corrections and hence all higher-order harmonics of the 
orbital frequency consistent with that PN order. The FWF is to be contrasted with the 
so-called restricted waveform (RWF) which retains only the leading-order harmonic at twice 
the orbital frequency. In applications to data analysis both the FWF and RWF should 



incorporate the orbital phase evolution up to the maximum aval 



which is currently 3.5PN 







14 



151 ] . Previous investigations [16 



able post-Newtonian order 



17 



FWF up to 2.5PN order 2 . Recently Kidder 



181 ] have obtained the 



pointed out that there is already enough 



information in the existing PN results [17| to control the dominant mode of the waveform, 
in a spin-weighted spherical harmonic decomposition, at the 3PN order. This mode, having 
(£,m) = (2,2), is the one which is computed in most numerical simulations, and which is 



1 As usual, we refer to nPN as the order equivalent to terms ~ (v/c) 2n in the asymptotic waveform (beyond 
the Einstein quadrupolc formula), where v denotes the binary's orbital velocity and c is the speed of light. 

2 The computation of the FWF is more demanding than that of the phase because it not only requires 
multipolc moments with higher multipolarity but also higher PN accuracy in many of these multipolc 
moments. This is why the FWF is known to a lower PN order than the phase. 



2 



therefore primarily needed for 



shall extend the works 



16 



17 



comparison with the PN waveforms. In the present paper we 



19j by computing all the spin-weighted spherical harmonic 



modes (£, m) consistent with the 3PN gravitational polarisations. 



The data analysis of ground-basec 



on the RWF approximation 



21 



and space-based detectors has traditionally been based 



22 



23 



24 



25 



261 ] . However, the need to consider 



the FWF as a more powerful template has been 



more accurate parameter estimation 



and the detection rate 



31 



27 



28 



29 



emphasized, not only for performing a 



32 



301 ]. but also for improving the mass reach 



331 ] . Another motivation for considering the FWF instead 



of the RWF is to perform cosmological measurements of the Hubble parameter and dark 
energy using supermassive inspiralling black-hole binaries which are known to constitute 
standard gravitational- wave candles (or sirens) in cosmology 34j, |35(. Indeed it has been 
shown that using the FWF in the data analysis of LISA will yield substantial improvements 



(with respect to the RWF) of the angu 
distance of gravitational- wave sirens [36J, 



ar, resolution and the estimation of the luminosity 
This means that LISA may be able to uniquely 



identify the galaxy cluster in which the supermassive black-hole coalescence took place, and 

the source which is crucially needed for 



thereby permit the measurement of the red-shift o 
investigating the equation of state of dark energy . 

It turns out that in order to control the FWF at the 3PN order we need to 
the multipolar post-Minkowskian (MPM) wave generation formalism 



38 



39, 



urther develop 



40 



41 



42 



43]. 



The MPM formalism describes the radiation field of any isolated post-Newtonian source and 
constitutes the basis of current PN calculations 3 . In this formalism, the radiation field is 
first of all parametrized by means of two sets of radiative multipole moments [47J. These 
moments are then related (by means of an algorithm for solving the non-linearities of the field 
equations) to the so-called canonical moments which constitute some useful intermediaries 
for describing the external field of the source. Finally, the canonical moments are expressed 
in terms of the operational source moments which are given by explicit integrals extending 



17 



4S 



49J most of the 



over the matter source and gravitational field. In previous studies [131 . 
required source moments in the case of compact binaries were computed, or techniques were 
developed to compute them. The important step which remains here is to refine, by applying 



An alternative formalism called DIRE has been developed by Will and collaborators 



44 



45 



3 



the MPM framework, the relationships between the radiative and canonical moments — this 
means taking into account more non-linear interactions between multipole moments — and 
between the canonical and source moments. The latter relationship involves controlling the 
coordinate transformation between two MPM algorithms respectively defined from the sets 
of canonical and source moments. 

The plan of this paper is as follows. In Section [TT] we recall the basic formulas for defining 
the FWF in terms of radiative multipole moments. Sections II I II and IIVI apply the MPM 
formalism to obtain general formulas for relating the radiative moments to the source mo- 
ments via the canonical moments. Section |V] summarises the results for all the relevant 
moments parametrizing the FWF at 3PN order. The time derivatives of source moments 
are investigated in Section |Vj and the various hereditary contributions are computed in 
Section IVHI The complete polarization waveforms at 3PN order are given in Section IVIIII 
for data analysis applications. Finally, the spin-weighted spherical harmonic modes of the 
3PN waveform are provided in Section [IX] for use in numerical relativity. 

For the benefits of readers we provide in Appendix El a list of symbols used in the paper 
together with their main meaning. 

II. THE POLARIZATION WAVEFORMS 

The full waveform (FWF) propagating in the asymptotic regions of an isolated source, 
hj^, is the transverse-traceless (TT) projection of the metric deviation at the leading-order 
1/R in the distance R — |X| to the source, in a radiative-type coordinate system = 
(cT, X). The FWF can be uniquely decomposed [47| into radiative multipole components 
parametrized by symmetric-trace-free (STF) mass-type moments Ul and current-type ones 
Vl. 4 The radiative moments are functions of the retarded time Tr = T — R/c in radiative 

4 The notation is: L = i\ ■ ■ ■ ig for a multi-index composed of I multipolar spatial indices i\, • • • ,iz (ranging 
from 1 to 3); similarly L — 1 = i\ ■ ■ -ii-i and aL — 2 = ai\ ■ ■ -ie—2] Nl — ■ ■ ■ Ni e is the product 
of £ spatial vectors Ni (similarly for Xl = ■ ■ -Xi e ); 8l = di 1 ■ ■ -di t and say d a h-i = d a di x ■ ■ -9i f „ 2 
denote the product of partial derivatives di — d/dx 1 ; in the case of summed-up (dummy) multi-indices 
L, we do not write the I summations from 1 to 3 over their indices; the STF projection is indicated 
using brackets, T(l) — STFfZj,]; thus Ul = U^l) and Vl — V(l) f° r STF moments; for instance we write 
x^iVj) — 7}(xiVj +XjV{) — j^SijX ■ v; e a h c is the Levi-Civita antisymmetric symbol such that £123 = 1; time 
derivatives are denoted with a superscript (n). 
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coordinates. By definition we have, up to any multipolar order 

+00 



TT 



AG 



^Z(N)J2^{Nl~2U ml „ 2 (T r ) --^—^N aL _ 2 e ab{k V l)bL ^ 2 (T R )j 



(2.1) 



Here N = X/i? = (AQ is the unit vector pointing from the source to the far away detector. 
The TT projection operator in ( 12. II) reads Vj^ = VikVji — \Pi{Pu where Vij = — NiNj 
is the projector orthogonal to the unit direction N. We introduce two unit polarisation 
vectors P and Q, orthogonal and transverse to the direction of propagation N (hence Vij = 
PiPj + QiQj)- Our convention for the choice of P and Q will be clarified in Section rVIIII 
Then the two "plus" and "cross" polarisation states of the FWF are defined by 

h 

1 1 ( T / - i . , 1 11 n 



c 2 R 



) 1 f 
Em{ Nl -* u ^ t ^- 



c(£ + iy 



(2.2) 



Although the multipole decompositions (12. ip and (12.21) are all what we need for our 
purpose, it will also be important, having in view the ongoing comparisons between the PN 



and numerical results 



3 impo 

By 



11 



121 ]. to consider separately the various modes (£, m) of the 



FWF as defined with respect to a basis of spin-weighted spherical harmonics. To this end 



we decompose h + and h x in the standard way as (see e.g. 

+00 i 



120) 



h + -ih x =J2 E h em Y^(&^) 



(2.3) 



1=2 m=-l 

where the spin-weighted spherical harmonics of weight —2 is function of the spherical angles 
(O, $) defining the direction of propagation N, 5 and is given by 



1 -2 



9f -\- 1 

^^4 m (9)e im * 

47T 



^ (-) fc + m)\(£ - m)\(£ + 2)\(£ - 2)! 



e 

cos — 

k\ (k - m + 2)\(£ + m - k)\(£ - k - 2)\\ 2 



2£+m-2k-2 



e 

sin — 

2 



(2.4a) 

2fc-m+2 

(2.4b) 



5 For the data analysis of compact binaries in Section I Villi the direction of propagation will be defined by 
the angles (0, $) = (i, ?) where i is the inclination angle of the orbit over the plane of the sky. 
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Here k% = max(0, m — 2) and hi = min(£ + m, £ — 2). Using the orthonormality properties 
of these harmonics we obtain the separate modes h im from the surface integral 



h 



h, — ihs 



3^(6,*) 



(2.5) 



where the bar or overline denotes the complex conjugate. On the other hand, we can also, 
following relate h lm directly to the multipole moments Ul and Vl- The result is 6 

G 



h 



tin 



u 



till 



-V 

c 



(2.6) 



where U £m and V are the radiative mass and current moments in standard (non-STF) 
□ 

guise [19J. These are related to the STF moments by 



i\ V 2£(£ - 1) 



-ylm ^_ 



+ 2: 



t\ ]l 2(£ + !)(£- 1) L L ' 



(2.7a) 



(2.7b) 



Here a^ 71 denotes the STF tensor connecting together the usual basis of spherical harmonics 
Y lm to the set of STF tensors iVm = N/^ ■ ■ ■ N i( \ (where the brackets indicate the STF 
projection). Indeed both Y lm and A^m are basis of an irreducible representation of weight 
£ of the rotation group. They are related by 

i 



(6,$) = al n Y tm {Q^) 



(Off _L 1 \\\ 

Y £m (Q, «&) = 1 J/" a? 1 N {L) (&, $) , 



(2.8a) 
(2.8b) 



with the STF tensorial coefficient being 7 



ai m = / dQN (L) Y 



■z-z-im 



(2.9) 



As observed in [19j] this is especially useful if some of the radiative moments are known to 
higher PN order than others. In this case the comparison with the numerical calculation for 
these individual modes can be made at higher PN accuracy. 



We have an overall si gn d ifference with due to a different choice for the polarization triad (N, P, Q). 



7 The notation used in 



19 



47| is related to ours by y£ 



m _ (21+1)!! —Ira 



47rf! 



Q 7 
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III. RELATION BETWEEN THE RADIATIVE AND CANONICAL MOMENTS 



The basis of our computation of t 



Minkowskian (MPM) formalism 



38 



39 



re radiative moments is the multipolar-post- 



40 



41 



42 



43j which iterates the general solu- 



tion of the Einstein field equations outside an isolated matter system in the form of a 
post-Minkowskian or non-linearity expansion. The formalism is then supplemented by a 
matching to the PN gravitational field valid in the near zone of the source. In this Section 
and the next one we sketch the main features of the MPM iteration of the exterior field 
while limiting ourselves to quadratic non-linear order because this is what we need for the 
new terms required in the FWF at 3PN order 8 . We shall work with harmonic coordinates 



(ct,x), which means that 



8^ = 0, (3.1) 

where the "gothic" metric deviation reads h al3 = sj—g g al3 — i] al3 , with g the determinant 
and g af3 the inverse of the usual covariant metric, and with r] a " = diag(— 1, 1, 1, 1) being an 
auxiliary Minkowskian metric 9 . Up to quadratic non-linear order the vacuum Einstein field 
equations take the form 

Uh aP = Nf(h) + 0(h 3 ) , (3.2) 

where □ = rfd^dv is the flat space-time d'Alembertian operator, and where N% denotes 
the quadratic part of the gravitational source term in harmonic coordinates — a quadratic 
functional of h and its first and second space-time derivatives given explicitly by 



Nf (h) 



h'^d^h^ + -8 a h^8 l3 h^ - ^d a hd p h 



2d {a h fM/ d fl h^ u + duff* {V h? + d^ u ) 



+ T] 



ai3 



-d p h^d P h^ + \d^hd P h + \d^h vp d U W P 

4 o 2 



(3-3) 



with h = jf^hfa,. The four-divergence of this source term reads 



-8 a hDh+ - 
4 2 



d a h 



(3.4) 



8 Cubic non-linearities do contribute at the 3PN order in the form of "tail-of-tails" but those have already 
been computed [5o| . 

9 Beware of the fact that the TT waveform defined by (|2.1[) differs by a sign from the spatial components 
of the gothic metric deviation, hj^ = —Vj^h^i + 0(h 2 ). 



In this paper we shall consider two explicit constructions of the quadratic-order external 
metric following the MPM formalism. The first construction, dealt with in this Section, 
will be parametrized by two (and only two) sets of moments, mass moments Ml and cur- 
rent moments Sl, which are referred to as the canonical multipole moments. The canonical 
moments are crucially distinct from the radiative moments Ul and Vl, and the MPM con- 
struction will provide the relations linking them to Ul, Vl- The second construction (in 
Section HV|) will deal with the link between Ml, Sl and six sets of moments II, Jl, Wl, Xl, 
Yl and Zl collectively named the source multipole moments. Among these, the moments II 
(mass-type) and Jl (current-type) play the most important role, while for reasons explained 
below the other moments Wl, Xl, Yl and Zl are called the gauge multipole moments and 
will appear to be subdominant. 

Armed with such definitions, the computation of the radiative field (j2.1j) - (l2.2p proceeds 



in a modular way (see Section 6 of [43] for further discussion). We start with relating the 
radiative moments {Ul, Vl} to the canonical moments {Ml,Sl} which are to be viewed 
as convenient intermediate constructs relating the radiation field and the matter source. 
The canonical moments are then in turn connected to the actual multipole moments of the 
source {II, Jl, Wl, Xl, Yl, Zl}. The point of the above strategy is that the source moments 
admit closed-form expressions as integrals over the stress-energy distribution of the matter 
and gravitational fields. The expressions of II, ■ ■ ■ ,Z L for general sources are given by 
(5.15)-(5.20) in 43( and shall not be reproduced here. 10 Note that the above formalism 
can be applied only to PN sources, which remain confined in their own near zone; the final 
expressions of the source moments are valid only for sources that are semi-relativistic like 
inspiralling compact binaries. 

Consider the so-called "canonical" construction of the MPM metric in harmonic coordi- 
nates, designated that way because it is based on Thome's 4?| canonical expression for the 
linearized approximation h^ n x [given by (I3.6P below] . The MPM metric is parametrized by 
the canonical multipole moments Ml and Sl and reads, to quadratic order, 

KL = Gh a J Ql [M L , S L ] + G 2 h a c L[M L , S l ] + 0(G 3 ) , (3.5) 



10 Below we give the source moments needed at the 3PN order in a form already reduced to the case of 
compact binaries in circular orbits. 

8 



where the Newton constant G is introduced as a convenient book-keeping parameter for 
labeling the successive non-linear approximations. From (I3.1l) - (l3.2p the linearized approx- 



imation ft-cani obviously satisfies d fl h^ nl = together with Dh 



a/3 
can 1 



0. Following 



38, 



39| 



we adopt the following explicit retarded solution of these equations, 

d L [r~ l M L (t - r/c) 



h 00 

' 4 can 1 



_1 V 



(3.6a) 



h 0i - — 

"-canl c 3 



it-i(* - r / c )J + j£iabdaL-i [r 1 S bL - 1 (t - r/c 



(3.6b) 



9, 



i-2 



r- X M^ L _ 2 {t-r/c) 



+ 



2i_ 

TV 



■d, 



aL-2 



r 1 £ab{iSfJ bL _ 2 {t-r/c) 



(3.6c) 



with F^ n \t) denoting n time derivatives of F(t). These expressions represent the most 
general solution of the vacuum linearized field equations modulo a change of gauge [471 ] . 



Next, the quadratically non-linear term h^ n2 — and subsequently all non-linear terms 



is constructed by the following algorithm. We first define 



u 



a/3 
can 2 



B=0 



r 



(3.7) 



where represents the ordinary (flat) retarded integral operator acting on the source 
N 2 (^-cani) which is obtained by insertion of the linearized metric (I3.6P into the quadratic 
source term given by (13.31) . The symbol FP#=o refers to a specific operation of taking 
the finite part when the complex number B tends to zero. Such a finite part involves the 
multiplication of the source term by a regularization factor (r/r ) B , where r represents an 
arbitrary constant length scale (and B G C). The finite part is necessary for dealing with 
multipolar expansions which are sin gula r at the origin r = [like in (13.61) ]. It will not be 



further detailed here and we refer to 



38 



41 



43| for full details. The point is that the object 



(13. 7p obeys the d'Alembertian equation we want to solve, namely 



0^2 = Nf(h canl/ 



(3.8) 



However, such a solution is a priori not divergenceless and so the harmonic coordinate 
condition needs not to be satisfied. To obtain a solution which is divergenceless we add to 



9 



/ M li can2' 



M can2 another piece v^ n2 defined as follows. Computing the divergence w" SLn2 = d u u^ ■■>■ «'c 
readily find 



W can2 — U rct 
r>=0 



r J r 



(3.9) 



where we used the fact that the source term of (I3.8p . as an immediate consequence of 
(13.4p . is divergenceless, dpN^ n2 = 0- Again, because the source term is divergenceless, the 
divergence u>" an2 must be a (retarded) solution of the source- free d'Alembertian equation, 
□u^ an2 = 0. This can also be checked from the fact that there is a factor B explicit in the 
source of (13.91) (appearing because of the differentiation of the regularization factor r B ), and 
therefore the finite part at B = is actually equal to the residue in the Laurent expansion 



38|. 



when B — > 0, and is necessarily a retarded solution of the source-free equation 

Given any vector of the type w" an2 , i.e. one which is of the form of a retarded solution 
of the d'Alembertian equation, we can always find four sets of STF tensors Nl, Pl, Ql and 
Rl such that the following decomposition holds, 



W can2 



(3.10a) 



w 



can 2 



Y,dL[r- l N L {t-r/c) 

1=0 

J2 d *L [r' l P L {t-r/c) 

1=0 

+oo 

+ {^-i V~ l QiL-i{t - r/c)] + e iah d aL ^ [r- l R hL ^{t - r/c)] } . (3.10b) 
t=i 

From this decomposition (which is unique) we define the object v^ n2 by the formulas 11 

-r-'N^ + d a [r- 1 (-N™ + - 3P«)] , (3.11a) 

(3.11b) 



..oo 

7 can2 



can 2 



V can2 



Q 



(-D + 3 p(D 



+ 00 s 
0—1 V 



+ d. 



'L-2 



r'\N^ L _ 2 + 3P®_ 2 - Q tjL - 2 )] - 1daL-2 [r^s ab{l R 3) bL-2] j • (3.11c) 

The superscript (— p) denotes the time anti-derivatives (ie. time integrals) of the moments. 
Such anti-derivatives yield some secular losses of mass and momenta by gravitational radia- 
tion which have been checked to agree with the corresponding gravitational radiation fluxes, 



ii 



We are adopting here a modified version of the MPM algorithm (with respect to [38|) as proposed in [51 1 
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see e.g. (4.12) in 5jj. The formulas (13.111) have been conceived in such a way that the 
divergence of the so defined v^f n2 cancels out the divergence of u^ n2 which is w^ an2 - I* 1 the 
following we shall denote by V a/3 the operation for going from a vector such as (I3.10p — a 
retarded solution of the source-free wave equation — to the tensor (13. lip . We therefore pose 

<fn 2 = V a >can 2 ] , (3.12) 

and as mentioned before this tensor immediately satisfies Dv^ n2 = (which is obvious) and 
also 

5 Xan2 = -<an2- ( 3 - 13 ) 

This property can be directly checked from (13.111) and (I3.10p . Finally, it is clear from (13. 8p 
and (I3.13P that by posing 

,a/3 a/3 , a/3 (3 14.) 

"-can 2 ~~ "can 2 ' y can2 > K - 1 ^) 

we solve the Einstein vacuum field equations at quadratic order, namely 

□^2 = ^2, (3.15a) 
<VC 2 = 0. (3.15b) 

The MPM algorithm can be extended to any post-Minkowskian order n. 

The structure of the quadratic metric h^ n 2 so constructed has been investigated in pre- 



vious works 



40 



51] . It consists of two types of terms: those which depend on the source 
moments at a single instant, namely the current retarded time t — r/c, referred to as in- 
stantaneous terms, and the other ones which are sensitive to the entire "past history" of the 
source, i.e. which depend on all previous times (r < t — r/c), and are referred to as the 
hereditary terms. The hereditary terms are themselves composed of three types of contribu- 
tions, the tail integrals — made from interaction between the mass of the source M and the 
time-varying moments and Sl (having i > 2) — , the memory integrals responsible for 



the so-called "non-linear me mory " or Christodoulou effect 52j, |53|, |54| (investigated within 



the present approach in 



40 



5jj), and semi-hereditary integrals which are in the form of 
simple anti-derivatives of instantaneous terms and are associated with the secular variations 
of the mass, linear momentum and angular momentum. The semi-hereditary integrals are 
given by the time ant i- derivatives present in the formula (13.111) . 

11 



To obtain the radiative moments we expand the metric at future null infinity in a radiative 
coordinate system = (cT,X l ), which is such that the metric admits an expansion in 
simple powers of 1/R without the logarithms which plague the harmonic coordinate system 



x^ = (ct,x l ) [39|. Up to quadratic order and for all multipole interactions we consider, 
we find that it is sufficient to define for the radiative coordinates X' 1 = x % and (denoting 
T„ = T- R/c) 

T a = t - r -- 2 -£E 1 jL) +0{ ^), (3.16) 



c c 3 \r 

where tq is the length scale introduced in ( 13. 71) . Expanding the metric when R — > oo 
with Tr = const, and applying the TT projection we obtain the radiative moments Ul 
and Vl we are seeking by comparing with their definition in (12. ip . At linear order the 
radiative moments agree with the £-th time derivatives of the canonical moments, and 
S L e \ At quadratic order we find that tail and non- linear memory terms appear; these have 



already been investigated in 



50 



5 11] . 12 Their general structure will also be given in (15 .ip 
below. Finally, we have numerous instantaneous terms whose determination necessitates the 
straightforward but long implementation of the MPM algorithm (I3.7p -( l3~14j) . This is the 
work required here: we have implemented the MPM algorithm in a Mathematica program to 
obtain all the instantaneous terms needed to control the 3PN waveform. The presentation 
of the results is postponed to Section IV A[ 

IV. RELATION BETWEEN THE CANONICAL AND SOURCE MOMENTS 
A. General method 

We next need to connect the canonical moments {Ml, Sl} to a convenient choice of mo- 
ments that are suitably defined to play the role of source moments. As it turns out, the 
source moments are best represented by six multipole moments {II, Jl, Wl, Xl, Yl, Zl} ad- 
mitting closed-form expressions in the form of integrals over the source and the gravitational 
field. To define them we consider a MPM construction which is more general than the one 



12 The semi-hereditary integrals associated with secular gravitational radiation losses do not contribute to 
the radiative moments. 
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given by (13.51) . namely (still up to quadratic order) 

KL = Gti£*i[Iu Jl, W l ,---, Z l ] + G 2 h a J n2 [I L , J L , W L , - • • , Z L ] + 0(G 3 ) , (4.1) 

where the linearized metric h^ n 1 is defined by the canonical expression h^ n 1 explicitly given 
in (I3.6p but parametrized by {II, Jl} instead of {M L , Sl}, and augmented by a linearized 
gauge transformation associated with some vector (pf parametrized by the remaining mo- 
ments {Wl,Xl, Yl, Zl} which can thus rightly be called the gauge moments. Thus, 

Ci = h cL[h, Jl] + d<p?\W L , X L , Y L , Z L ] , (4.2) 

where for any vector ip® we denote the gauge transformation by 

dtp? = d a p{ + dPtf - r! a %^ . (4.3) 

Note that d^dip? = Op". The expression of ipf in terms of the gauge moments is 

& E [r- l W L (t - r/c)} , (4.4a) 

d = ~ ? E {j w d ^ V' lx ^ - r / c )] ( 4 - 4b ) 

~ i E { dL - 1 V~ XY ^ ~ r / c )] + -Ar^daL-x [r- l Z bL - X {t - r/c)] X . 

(4.4c) 

The quadratic metric h^ n2 will now be defined by the same algorithm as for the canonical 
metric in Section [Till but starting from the general linearized metric (14.2ft . The result will be 
another MPM metric (both the canonical a nd g eneral metrics are legitimate to describe the 



exterior field of any isolated matter source 



38] ) and we shall look for the relation between 



{M L , Sl} and {II, Jl, Wl, Xl, Yl, Zl} which is necessary in order that these two metrics 
differ by a coordinate transformation (at quadratic order), and therefore describe the same 
physical matter source. 

To proceed, we have to define 



u 



a/3 

gen 2 — '—'ret 



fp nzl 



-Y Nf(h gcnl ) 



(4.5) 



The only difference with ( 13.71) is that the quadratic source (I3.3P is computed from h^ nl 
instead of h^ nl . However, since the two linear metrics h^ nl and h^ nl differ by the gauge 

13 



transformation (14. 3[) the difference between the corresponding sources must have a specific 
structure, and we find 



Nf (h genl ) = N^(h c&nl ) + nnr + dA a 2 



r a/3 



a/3 



a/3 



(4.6) 



We employ the notation f!4.3[) for the gauge term dA^ 13 . The expressions of the tensor Q 
and vector A?? are determined with the help of (13.31) and read 



>a/3 



a/3 



can 1 



^2 



-/Ci»Avi+^(^i) 

As a consequence of (14. 7p we easily verify that 



(4.7a) 
(4.7b) 

, M ** 2 ^ 2 (4.8) 

This relation is consistent with the fact that the source term N 2 ^ is divergenceless [because 
of (13.41) ]. Hence we see that the divergence of (14.61) is automatically verified, where we use 



W + A, 



0. 



the fact that d^dA 2 ^ 



□Af. 



Applying our specific finite part of the retarded integral operator on both sides of H4.6[) 



we obtain the relation between u^ n2 defined by (14.51) and the corresponding u^ n2 defined 
by (13. 7p in the canonical algorithm, namely 



u 



a/3 



a/3 ^a/3 
"can 2 ^ i4 2 



''gen 2 — "can 2 T "2 T ^^2 T ^2 1 2 ■ ( 4 -9) 

The difference between the two prescriptions is made of various terms. The terms Q 2 ^ 



and 



represent what we would expect if the operation of taking the finite part of the 



retarded integral would commute with partial derivatives. Here the gauge transformation is 
associated with the gauge vector defined by the finite part of the retarded integral of A 2 , 



FP □ 

B=0 



-1 



B 



A a 

^2 



(4.10) 



The last two terms X 2 ^ and Y 2 come from the non commutation of the finite part of the 
retarded integral operator FPs=o □~t(r/Yo) B with the differential operators □ and d which 
are present in front of the last two terms of (13.31) . respectively. We have 



X. 



a/3 



fp n~l 

B=0 rct 



2 i 



(4.11a) 



14 



v a(3 
I 2 



B=0 



dAf 



which can also be seen more formally as the action of "commutators" namely 



X. 



all 



Y, 



a/3 



FPD-t 1 , □ 

FPD 7eL 9 



a/3 
2 > 



^2 



(4.11b) 

(4.12a) 
(4.12b) 



Our notation for the commutators involved and for the partial derivative d should be 
clear. Here FP is a short hand for FPb =0 n~ e \(r/r ) B and we have used the fact that 
□ (FP n~]f) = f. It is evident that the non commutation of FP n ret with partial derivatives 
comes from the presence of the regularization factor r B . Thus X 2 ^ and Y 2 a/3 are built from 
the spatial differentiation of r B , i.e. dir B = Bn i r B ~ 1 , and will involve an explicit factor B 
in their sources. Their expressions read as 

, 2 3 = FT D~l B I ' ' 1 



Y 



af3 



B=0 



FP □ 

B=0 



B + 1 p 2 onf 

o ^2 o 

r or 



r- 



-i 



B (L\ 3 V± (-S^A* - 5 i/3 A« + r) al3 Ai 



(4.13a) 



(4.13b) 



In Section IIVBI we shall present a practical method to evaluate X 2 an d Y 2 13 at the lowest 



ra/3 



PN order, given the general quadratic-type structure for the source terms (j4.7aj) and (14.7bl) . 

The first part of the MPM algorithm u^ n2 has been obtained in (14. 9p . and we look now 
for the second part v^ n2 . To this end we compute the divergence u>g en2 = °^u^ n2 . Using 
(14. 9 ft and the property ( 14. 8ft we readily find that 



w 



gen 2 



W. 



can 2 



(4.14) 



where we pose for simplicity 

Uf = Xf + Yf . (4.15) 

The structure (14.131) of X 2 ^ and Y 2 a/3 involving the retarded integral of a source term contain- 
ing an explicit factor B implies that U 2 is necessarily a retarded solution of the source-free 
d'Alembertian equation, UU 2 = 0. Hence, there must exist ten STF tensors A L ,B L , - ■ ■ ,L L 
(functions of the retarded time) parametrizing the ten components of U 2 in such a way 
that 



(4.16a) 
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+oo 



U^ = Y,^l [r~ l B L {t-r/c 

£=0 

+00 

+ E i^- 1 [ r_lc «-i(* " r /°)} + ^abdai-i [r-'D^t - r/c)] } , (4.16b) 
+00 

U 2 = E - r / c )] + *A [ r ~ lF L(t - r/c)] } 

+00 

+ E [ r_1 Cj)L-i(t - r/c)] + E ab{ id j)aL -i [r^HtL-iit - r/c)] } 

e=i 
+00 

+ i dL ~ 2 V^ K ^-2it - r/c)} + <9 ai _ 2 [r- 1 £ ab(i L J>i _ 2 (t - r/c)} } . (4.16c) 



The divergence of this tensor, W% = d^U 2 , will also be of that form and hence there will 



exist four STF tensors N' L , ■ • • , R' L such that 



W " = J2 dL V' 1 ^ - r / c )] > ( 4 - 17a ) 
+00 

W* = J2 d ^ [r- l P' L (t-r/c)} 

+00 

+ E {^-1 ['•"ViL-iC* - r/c)] + e iab d aL -t [r-'R'^t - r/c)} } . (4.17b) 



1=1 



The four tensors N' L , ■ ■ • , i?^ play exactly the same role as Nl, • ■ ■ , Rl in (13.101) . and we 
shall apply the same algorithm as the one going from (I3.10p to (13. lip . Thus, we define from 
the components of W£ a new tensor V% by this algorithm, which was denoted by V a/3 in 
f[3TT2|) . Hence 



so that in component form this tensor reads 



VT = -r-'N'^ + d a [r- 1 (-Nt 1] + Q'^ ~ 3^)] 



K 0i = r 



1 -Q?- 1 ) + 3if> )-e iab d a r-iR'^l - J>_x l^L-i] 



+00 



(4.19a) 
(4.19b) 



+ OO ,- 



+ 3 



L-2 



r _1 (iV{L + 3^}2_ 2 - QU- 2 ) - 2d aL - 2 [r-^i?-)^] }> . (4.19c) 



y(2) 
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However, in the present case the tensors N' L , • • • , R' L can be directly related to the ones 
parametrizing (14.161) . By computing the divergence W£ = d^U^ we readily find 

N' L = A^ + Bf ] + C L) (4.20a) 
P' L = Ef + F L + l -G L + 4 1} , (4.20b) 

Q'l = \Gf +K L + cf , (4.20c) 

R'l = \H? + \L L + Z>« • (4.20d) 

Thus V 2 °^ can be expressed directly in terms of Al, ■ ■ ■ , Ll by substituting ( 14.201) into (14.191) . 
In doing so we shall discover that the time anti-derivatives present in ( I4.19P become in fact 
"instantaneous" because they are cancelled by some time derivatives coming from (14.201) . 
By construction of (14. 19j) we have at once UVf = and dyY^ = -W£. Applying the 
MPM algorithm we therefore find for the second part of the algorithm, 

v a J^ = v a c L 2 + Vf. (4.21) 

Gathering the results (14.91) and (I4.2ip the complete quadratic-order metric is obtained as 

KL2 = KL + + d<j>f + Uf + V? , (4.22) 

and it satisfies the vacuum Einstein field equations in harmonic coordinates, i.e. 

□/£!■, = JV£ 2 , (4.23a) 
<Vg n2 = 0. (4.23b) 

To find the relation between the source and canonical moments we notice that the sum of 
the last two terms in (I4.22[) is a solution of the linearized vacuum equations, since it satisfies 
□ (C/2 + ^2 ) = and also d^U^ + V^) = 0. It must therefore be of the form of the 
general solution h^ nl of these equations which has been given in ( 14. 21) . i.e. there should 
exist some moments 5Il and 5Jl representing specific corrections to II and (necessarily 
at quadratic order) and some gauge vector ijj^ sucn that 

ra/3 , rra/3 _ h a(3 r r r r t i . o , a/3 



Ur + = hZ 1 [5I L , SJ L ] + tyf . (4.24) 
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Let us prove that the corrections we seek to the moments II and Jl that are needed to 
reproduce the canonical moments are indeed provided by these SIl and SJl, i.e. 

M L = I L + G SI L + 0{G 2 ) , (4.25a) 
S L = Jl + G SJ L + 0(G 2 ) . (4.25b) 

To this end we have to check that the general metric h^JI^, Jl, Wl, ■ • • , Z L \ constructed 
at quadratic order in (14.11) is isometric — differs by a coordinate transformation — to the 
canonical metric h^ n [ML, Sl] given by (13.51) . This immediately follows from (I4.22p and 
( I4.24p which permits us to recast the general metric (14. ip into the form 

a @ r /l/f c 1 _L f) a @ 



h£[I L ,J Lt .--] = G Ci[^,^] + ^ 



+ G 



h a J n2 [M L , S l ] + Qf + dtp? + 0{G A ) , (4.26) 



i3 ^ 



where we have posed <^2 = 02 + ^21 an d where higher-order powers of G are consistently 
neglected. From this result we conclude that h^ n [I L , J L , ■ ■ ■} and h^ n [M L , Sl] differ by the 
coordinate transformation 

<en = <n + Gtf + G\« + 0(G 3 ) , (4.27) 

as we have recognized that represents precisely the quadratic non-linear part of that 
coordinate transformation, i.e. the term which makes it to differ from a linearized gauge 
transformation. Hence we have proved that the two sets of moments {II, Jl, Wl, X l , Y l , Z l } 
and {M L , Sl} related by (14.251) are physically equivalent — they describe the same physical 
matter source. Note that the relations (I4.25P give the canonical moments as functionals of 
the full set of source moments {II, Jl,Wl, Xl,Yl, Zl}. Consequently, just two moments 



Ml and Sl are still sufficient to describe the external field of any source 



381 ] . Notice also 



that M L and Sl are almost equal to II and Jl in the sense that the corrections Sir and SJl 



in (14.251) will turn out to be very small in a PN expansion, being of order 2.5PN [42j]. This 
is of course the result of the fact that the gauge moments {Wl, X l , Y l , Z l } do not play any 
physical role at the linear approximation, where the coordinate transformation reduces to 
the gauge transformation. However, since the theory is covariant with respect to non-linear 
diffeomorphisms and not merely with respect to linear gauge transformations, the moments 
{Wl, X l , Y l , Z l } do play a role at the non-linear level. 

18 



B. Practical implementation 



Finally let us sketch our practical method to compute the correction terms 8Il and 
5Jl. We remark first that they come via (14.241) from the ten STF tensors A^,--- ,Ll 
parametrizing as given by (14.161) . We can therefore express 51 l and 5Jl directly in 
terms of Al, • • • , L L by following in details the steps (I4.17p - (l4.20p . The result is 

< Yd. 



SI, = -c 2 



4 



A L + 4B { P + 3E { ? + 3F L + G L 



5Jr 



,(-)^ + l)! 



-H 



(i) 



(4.28a) 
(4.28b) 



The next problem is to compute the tensors A^, ■ ■ ■ , in the PN approximation. These are 
defined from the two objects X% and Y 2 which are given in particular by their commutator 
form (I4.12p . We thus need to compute the commutator between the operator FP and 
derivative operators, when applied either on the terms fij or ^2^- The relevant point for 
our purpose is that the general structure of these terms at the quadratic order is known. 
Namely fi^ and are made of quadratic products of retarded multipolar waves, i. e. are 
given by sums of terms of the type 



/Cpq = d {P) [r~ l F(t - r/c)] d {Q) [ r - l G{t - r/c)] 



(4.29) 



where the functions F and G stand for some time derivatives of moments in the list 
{I Li J L-, Wl, Xl, Yl, Zl}- It is convenient to suppress the indices on these moments and 
to write only the "active" indices appearing in the spatial multi-derivatives dp and 8q, com- 
posed with the multi-indices P = ai ■ ■ ■ a p and Q = bi ■ ■ ■ b g (p and q being the number of 
partial derivatives in dp and 3q). Furthermore the multi-derivatives in (14.291) are chosen 
to be STF (this can always be assumed modulo a possible STF decomposition), hence the 
brackets () surrounding their indices. The problem is therefore reduced to that of evaluating, 
in the PN approximation, the quantities 13 



y l p Q 



FPn-J, □ 



FPD-J,^ 



K-PQ 1 
K-PQ ■ 



(4.30a) 
(4.30b) 



13 In the case of y l pq we can restrict ourselves to a spatial derivative d z because the time derivative dt 
commutes with the operator FPD~J, thus 3*pg = 0. 
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Indeed and are given by some sums of terms of the type Xpq and ypQ respectively 
(and multiplied by appropriate constant tensors involving Kronecker symbols to perform the 
needed contractions). 

The term Xpq has in fact already been computed at the lowest PN order in the Appendix 
of 42|. The result turned out to be quite simple, namely 



X, 



PQ 



-d 



(PQ) 



8 P)0 FWG + 6 . a FGV> 



0,g 



o 



(4.31) 



Here the functions are evaluated at retarded time t — r/c (with F^ and denoting the 
time derivatives), and 5 Pi o and <$o, 9 denote the usual Kronecker symbols. As we see in (I4.3ip 
the two STF multi-derivative operators dm\ and d/Q) originally present in (14.291) have merged 
into a single STF derivative operator <9(pq) withp+g indices. The formula (14.311) constitutes 
a useful practical lemma for doing computations at the lowest PN order. Because of the 
factor 1/c in front of (14.311) the PN "parity" of the result (14.31 j) will be opposite to that of 
the source term (I4.29p . which in practice will typically be even. As a consequence we shall 
find that the PN order of is dominantly "odd", starting in fact with 2.5PN. 

As for the term ypq, it was not required in (42)] but will play a role here for the waveform 
at the 3PN order. We have worked out the equivalent of (I4.3ip for this term, and find, still 
at the lowest PN order, 



y 



PQ 



c(2p P + + 2 9 q + l){ 6p '° 5 ^ dpQ ^ 



+0 



(4.32) 



Consistent with our notation we write P — 1 = a± ■ ■ ■ a p _i and Q— 1 = b% • • • Again there 
is a factor 1/c and we shall find that the corresponding is dominantly "odd" , starting 
at 2.5PN order. Note that the new lemma ( 14.321) is not independent from the previous one 
(14.3ip and is actually more general than it. Indeed, by computing the divergence of y P Q 
using its definition (14.30bl) . we get 



FPD-J, □ 



K 



PQ 



FPDri, d l 



i'^PQ , 



(4.33) 



which can be used to check the consistency of the two formulas ( 14.311) and (14.321) . The 
results needed at 3PN order for the relation between the canonical and source moments as 
obtained by these means — namely the formulas ( 14.281) and the lemmas ( I4.31l) -( )4.32l) — are 
reported in Section IV Bl 
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V. THE MOMENTS FOR 3PN WAVEFORM 

Using the MPM algorithm of Section II III the radiative moments {Ul, Vl\ are related 
to the canonical moments {Ml, Sl}, and following Section HVl the canonical moments are 
in turn expressed in terms of the source moments {II, Jl, Wl, Xl, Yl, Zl}. In the current 
Section we present the results (skipping some details) of the computation of all the moments 
needed for controlling the FWF in the case of compact binary systems up to 3PN order. 

A. The radiative moments for 3PN polarisations 

To obtain the gravitational polarisations at 3PN order one must compute: the mass 
radiative quadrupole with 3PN accuracy; the current radiative quadrupole Vy and mass 
radiative octupole U^ with 2.5PN accuracy; mass hexadecapole U^i and current octupole 
V ijk with 2PN precision; U ijMm and V ijM up to 1.5PN order; U ijMmn , V mm at 1PN; U ijk i mno , 
Vijkimn at 0.5PN; and finally Uijki mn o P , Vijkimno to Newtonian order. The relations connecting 
Ul and Vl to the canonical moments Ml and Sl are first obtained following the MPM 
method of Section HTT1 14 

The quadratic contributions to the radiative mass (resp. current) moments are found in 
the form of sums of terms 5 2 Ul(u) (resp. S 2 Vl(u)/c) whose general structure reads 



The power of 1/c in front is chosen in such a way that m represents the PN order of our 
calculation of the waveform, i.e. m = 6 at the 3PN order. The capital letter A stands either 
for M or S, meaning that we are considering in (15.11) interactions between canonical moments 
of the type M^ l} M^f or M^ l} S^f or S^f , with the superscript (p) denoting time 
derivatives, and the multi-indices K\ and K 2 having length k\ and k 2 {e.g. K\ = a\ ■ ■ ■ a^). 15 
The kernel Xlk 1 k 2 has itself an algebraic structure made of a sum of products of Kronecker 
and Levi-Civita symbols. Its physical dimension depends on time only, and each of its 

14 We have implemented the MPM algorithm on the algebraic computing software Mathematica using the 

powerful tensor package xTensor [551 ]. 

15 The reasonning we shall make can be easily generalized to n-th non-linear terms 6 n UL, 5uVl/c involving 
n canonical moments Ag?\ • • • , A^ n \ 





c 



(5.1) 
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three sets of indices, L, K\ and K 2 , is symmetric and trace-free (STF). For instantaneous 
terms, which are functions of the multipole moments A^\ evaluated at the instant of 
emission u = t — r/c, it is proportional to the Dirac function S(s — u). 

The above structure does not exist generically for an arbitrary pair of multipole moments 
nor for any arbitrary value of k\ and k 2 . A closer look will actually allow us to reduce the 
number of terms in the source we shall focus on to a few ones, making our task much easier. 
As a product EijkE a b c can always be transformed into a linear combination of bia'Sjb'Skc' with 
{a',b',c'} = {a,b,c}, the number e of Levi-Civita symbols in each of the individual terms 
55...5ee...e composing Xlk 1 k 2 ma y De reduced to or 1. The symmetry of parity implies 
that this number is the same for all terms. Now, if e = 0, the integer £+ k\ + k 2 is even (equal 
to twice the number of Kronecker symbols) and all indices in L must contract with an index 
of K\ or K 2 . Thus, we must necessarily have ki + k 2 > t. On the other hand, if e = 1, the 
Levi-Civita symbol carries one index from each of the three STF sets, so that there remain 
I — 1 free indices of type L carried by Kronecker symbols, as well as k\ — 1 indices (resp. 
k 2 — 1) of type K\ (resp. K 2 ) involved in the contraction of some 5's with the multipole 
moments. Then, the same arguments as before show that (I — 1) + ik\ — 1) + ik 2 — 1) must 
be even with {k\ — 1) + {k 2 — 1) > I — 1. The previous constraints can all be summarized 
by the single statement that k\ + k 2 — i — e is always an even positive integer. 

The structure of the quadratic interactions may be further refined by noticing that only 
the multipole moments that have dimensions compatible with (15. ip are allowed to enter 
5 2 U, 5 2 V/c. Let us pose for later convenience [A^] = [M] [L]«i+*i-pi [V] a ^ and [ApJ] = 
[M] [L] a2+fc2 - p2 [V] a2+P2 , where [M], [L] and [V] denote the dimension of a mass, a length 
and a velocity. Equating [U L ] = [V L /c] and [dsxiK^] [G / c m ~ e+2 A^ A^f] on the one 
hand, remembering on the other hand that [xLKtK 2 ] is a certain power q G Z of the time 
dimension [T], we find: 

(dj + ati + ki) = m — 1 , (5.2a) 

8=1,2 

^(cti+pi) = q + m + 1. (5.2b) 

8=1,2 

Now, we know that k\ + k 2 — i — e 6 N. Moreover, the number e = 0, 1 of Levi-Civita 
symbols is itself governed by the parity symmetry. More precisely, defining the integers &i, 
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Q.2 and e, associated with A^J and 5 2 Ul, 5 2 Vl/c respectively, to be equal to zero when 

the latter multipole moments are of mass type or to 1 when they are of current type, the 
consistency of the transformation of both sides of equation ( 15. ip under parity imposes that 
e = dii + ^2 + e [mod 2]. As a result, the maximum multipolar order ki + k 2 — e of the 
radiative moments containing a quadratic interaction ^^jf is given by 16 



1 - ^ (of + oti 

4 = 1,2 



remainder 



on + e 



i=l,2 



(5.3) 



At the ^PN approximation, such a contribution exists only if 



> 2, or equivalently 



i=l,2V 



+ «j) + remainder [|(^ i=1 2 «j + e)] < m — 3. Once this necessary condition is 
fulfilled, the orders of multipolarity possibly affected by the piece of non-linear correction 
(pTT]) are 4ua(a», «», «i, e), £ max (ai, «», a i; e) - 2, • • ■ , 2 or 3, with e = (resp. e = 1) for 
mass (resp. current) radiative moments. The latter "selection" rules may be generalized 
to interactions of any post-Minkowskian order n, in which case m — 1 must be replaced by 
m+5— 3n in the expression (15.31) of £ max while all summation ranges become 1 < i < n. With 
the selection rules (I5.2p - (I5.3I) we are able to know beforehand which non-linear multipole 



interaction is needed to be computed in the radiative moments Ul, Vl at a given PIN 



3 



oraer. 



The result concerning the 3PN mass quadrupole moment Uij is already known [4d. 1501. 151| 
and we simply report it here. Actually, at 3PN order involves a cubically non-linear term, 
composed of the so-called tails of tails, whose computation necessitates an extension of the 
MPM algorithm to cubic order G 3 [50j. We have 



U ij {T R ) = MX>(T R ) 



(2),..- . . 2GM 



Tr 



dr 



In 



T, 



R 



T 



2r 



11 
+ 12 



G 



2 r Ta 
7 



a(i 



+-M^.M j)a - h -M^.M^ - -M^M?) + l -e ah{i MfS b \ 

j a{t J/ a j a(i j)a j a(i j/a g au \ l j)a " I 



n ( GM\ 

+2 (-) 



+ O 



dr 



In 2 



T, 



R 



2r 



70 V 2r 



+ 



124627' 
44100 



(5.4) 
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The remainder function means the usual division remainder: remainder [-j] = or 1 depending on whether 
N is an even or odd integer. 
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Notice the tail integral at 1.5PN order, the tail-of-tail integral at 3PN order, and the non- 
linear memory integral at 2.5PN. In the tail and tail-of-tail integrals, M represents the mass 
monopole moment or total mass of the binary system. The constant r in the tail integrals 
is given by tq = tq/c, where r$ is the arbitrary length scale originally introduced in the 
MPM formalism through ( 13. 71) . and appearing also in the relation between the radiative and 
harmonic coordinates as given by (I3.16p . 

The moments required at 2.5PN order are new with this paper (apart from the tails) 
and involve some interactions between the mass quadrupole moment and the mass octupole 
or current quadrupole moments. Which type of interactions is determined by using the 
selection rules discussed above. These moments are given by 17 



U ijk (T R ) = M$(T R ) + ^ ! TR dr 



In 



T r -t \ 97' 
2r J + 60 



+ 



G 



- -M,.M 



dr 



(3) „ #(3) 



l M§{T)M§ a {T) - -e ab{l M ( S{r)S^ b {r) 
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-M { fM { l\ +~M y y.M 

i rt 1 1 i £•* \ n 1 a nil 
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At 2PN order we have the standard tails and some previously known interactions of the 



17 In all formulas below the STF projection () applies only to the "free" indices denoted ijkl ■ ■ ■ carried by 
the moments themselves. Thus the dummy indices such as abc ■ ■ ■ are excluded from the STF projection. 
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mass quadrupole with itself 
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At 1.5PN we again have some non-linear interactions (new with this paper) involving the 
mass octupole and current quadrupole and given by 
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For all the other moments that are required, it is sufficient to assume the agreement between 
the radiative and canonical moments, 



u l {t r ) = m ( 1\t r ) + o 



v l {t r ) = s1\t r ) + o[^- 



(5.8a) 
(5.8b) 
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B. The canonical moments for 3PN polarisations 



Following the investigation of Section IIVI we now give the canonical moments in terms of 
source-rooted multipole moments. It turns out that the difference between these two types of 
moments — which is due to the presence of the gauge moments defined by (14.41) — arises only 
at the small 2.5PN order. The consequence is that we have to worry about this difference 
only for the 3PN canonical mass quadrupole moment M^-, the 2.5PN mass octopole moment 
Mijk, and the 2.5PN current quadrupole moment Sy. For the mass quadrupole moment, 
the requisite correction has already been used in and is given by 18 



M a = T a + — 



W^I tJ - jW] + o (±) , (5.9) 



where 1^ denotes the source mass quadrupole, and where W is the monopole corresponding 
to the gauge moments Wl (i.e. W is the moment having i = 0). At the PN order we are 
working, W is needed only at Newtonian order and will be provided in Section IV CI Notice 
that the remainder in (I5.9P is at order 3.5PN - - consistently with the accuracy we aim 
here. The expression (15. 9p is valid in a mass-centred frame defined by the vanishing of the 
mass dipole moment: J, = 0. Note that a formula generalizing (15.91) to all PN orders (and 
all multipole interactions) is not possible at present and needs to be investigated anew for 
specific cases. Thus it is convenient in the present approach to use systematically the source 
moments {II, Jl, Wl, Xl, Yl, Zl} as the fundamental variables describing the source. 

Similarly, the other moments and Sij will admit some correction terms starting 

at the 2.5PN order. We have computed these new corrections, together with recomputed 
and confirmed those in (15.91) . by following the method of Section IIVI i.e. evaluating the 
STF tensors Al, ■ ■ ■ ,Ll in (14.161) by means of the two practical lemmas (I4.3ip and (14.321) . 
then plugging these tensors into (14.281) . We also performed an independent calculation by 
implementing the general MPM algorithm of Section II III starting directly with the general 
linearized metric (14.21) parametrized by the source moments, instead of the canonical metric 
(13. 6p parametrized by the canonical moments. 

In this second approach, which fully confirmed the previous results, we need to know 



IS 



The equation (11.7a) in 13| contains a sign error with respect to the original result 42| (with no conse 



quence for any of the results in The correct sign is reproduced here. 
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beforehand the relevant multipole interactions and we used the same selection rules (I5.2p ~ 
(I5.3P as before. The only difference is that the quadratic interaction we consider, xA^\ 
are between any two source multipole moments composed of the main moments {II, Jl} and 
the gauge moments {Wl, Xl, Yl, Zl}, i.e. the letter A symbolizes now any of the /, J, W, 
X, Y or Z. By applying the selection rules ( 15.2j) -( pT3l) at 3PN order, i.e. for m = 6, it is 
straightforward to check that (i) no gauge multipole moment can enter cubic interactions up 
to our approximation level, and (ii) all quadratic interactions involving at least one gauge 
moment have to be instantaneous, meaning that q = — 1 for them. We can in fact determine 
all possible contributions by inspection. Their full list is given in table [H 



TABLE I: Non-linear corrections in and Sij involving at least one gauge multipole moment 
at 3PN order. The first entry indicates, for each interaction, which radiative moment it belongs 
to, whereas the second entry tells us how many time derivatives are involved. STF symbols are 
omitted. 
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in Mijk 




I^xWW IxWg> I$-*>xWj*> 




in 
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£iabl a j XW b 



Further rules of selection might be used to discard some candidates, but all the contri- 
butions to Uijk and Vij that are presented here have been computed explicitly. By retaining 
only the interactions that involve the pairs of multipole moments composing the elements 
of table [U 19 the source used in our algorithms could indeed be reduced to a finite, reason- 
ably small number of terms. However the detailed calculation of some of these interactions 

(2) 

turns out to yield zero; this is the case for instance of the interaction IxZh which does not 



contribute to S^. 

Finally our explicit results for and SV, are 
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(5.10b) 
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I.e. the interactions IxYijk, IiXYj, IxWijk, hjXWk, JiXYj, IxZy, JixWj and JijXW. 
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where Wi and Yi are the dipole moments corresponding to the moments Wl and Yl- The 
remainders in (15.101) are consistent with our approximation 3PN for the FWF. Besides the 
mass quadrupole moment (I5.9p . and mass octopole and current quadrupole moments (15.101) . 
we can state that, with the required 3PN precision, all the other moments Ml agree with 
their corresponding Jl, and similarly the Sl agree with Jl, namely 

1 



M L = I L + O 



s, Ji h CM I 



(5.11a) 
(5.11b) 



C. The source moments for 3PN polarisations 

We have finally succeeded in parametrizing the FWF entirely in terms of the source 
moments {II, Jl, Wl, X l , Y l , Z l } up to 3PN order. The interest of this construction lies 
in the fact that the source moments are known for general PN matter systems. They were 



obtained by matching t 
the source's near zone 



re external MPM field of the source to the internal PN field valid in 

Deen worked out in the case 
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43J. The source moments have 



of compact binary systems with increasing PN precision [13|, [17|, |48|, |49j . Here we list all the 
required J^'s and Jl's (and also the few needed gauge moments) for non-spinning compact 
objects and for circular orbits. We do not enter the details because the derivation of these 



13 



moments follows exactly the same techniques as in 

The only moment needed at the 3PN order is the mass quadrupole moment Lj, first 
computed for circular orbits in [13j and subsequently extended to general orbits in 49j . We 
write it as 

^,3 



vm I Ax(ij) + B 



Gm 



C 



Gm 



(i v 3) 



(5.12) 



The relative position and velocity of the two bodies in harmonic coordinates are denoted 
by x % = y\ — y\ and v % = dx 1 jdt = v\ — v\ (spatial indices are lowered and raised with the 
Kronecker metric so that X{ = x l and Vi = v l ). The distance between the two particles 
in harmonic coordinates is denoted r = |x|. The two masses are mi and m2, the total 
mass is m = m\ + m 2 (not to be confused with the mass monopole moment M), the 
symmetric mass ratio v = mxm^jm? satisfies < v < 1/4, and the mass difference ratio 
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is A = (mi — m 2 )/m which reads also A = ±\/l — 4z/ (according to the sign of m x — 777.2). 
To express the coefficients A, B and C in (15.121) as PN series we introduce the small post- 
Newtonian parameter 

7=^- (5.13) 



rcr 



With these notations we have (in the frame of the 'center-of-mass' and for circular orbits) 
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(5.14a) 



(5.14b) 
(5.14c) 



The coefficients A and i? correspond to conservative PN orders (which are even), while the 
coefficient C involves a single term at the odd 2.5PN order due to radiation reaction. 

Notice the appearance of logarithms in both A and B at the 3PN order. These loga- 
rithms have two distinct origins, depending on whether they are scaled with the constant 
r associated with the finite part prescription in (13.71) . or with an alternative constant de- 
noted r'o- The logarithms with ro will combine later with other contributions due to tails 
and tails-of-tails, and the constant ro will be absorbed into some unobservable shift of the 
binary's orbital phase, as can already be seen from the fact that ro is associated with the 
difference of origin of time between harmonic and radiative coordinates, see (13.161) . 

The other constant r' Q is defined by mlnr'o = mi Inr'i + m 2 lnr' 2 , where r'% and r' 2 are 
two regularization constants appearing in a Hadamard self-field regularization scheme for 
the 3PN equations of motion of point masses in harmonic coordinates The constant 

r'o is therefore present in the 3PN equations of motion and we shall thus also meet this 
constant in the 3PN orbital frequency given by (16 .4p below. The regularization constant r' 
is unobservable, since it can be removed by a coordinate transformation at 3PN order — r' 
can rightly be called a gauge constant. In practice this means that r'o will cancel out when 
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using the 3PN equations of motion to compute the time derivatives of the 3PN quadrupole 
moment, as will be explicitly verified in Section IVT1 20 

The list of required moments continues with the 2.5PN order at which we need the mass 
octupole and current quadrupole given by (with A = mi ~ m2 ) 
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Note also that the 3PN quadrupole moment 13|, |49( depended originally on three constants £, k, £ (called 



ambiguity parameters) reflecting some incompleteness of the Hadamard self-field regularization. These 



constants have been computed by means of the powerful dimensional regularization 15|, |58| , and we have 



replaced the result, which was £ 



9871 
9240 ' 



and C = -X, back into (|5. 14[) . 
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At 1.5PN order: 
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At 0.5PN order: 
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At Newtonian order: 
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(5.18b) 



(5.19a) 
(5.19b) 

(5.20a) 
(5.20b) 



The 2.5PN correction terms in 1^ and Jy, the 2PN terms in and Jij k , and the 1PN 
terms in lij k \ m and are new with this paper. The higher-order Newtonian moments 
1 ijklmno an d Jijkimn were also not needed before, but Newtonian moments are trivial and are 
given for general £ by 



I L = um s e (v) x {L) + O (J^j , 

J L _ X = Vm Sl {y) eab{i t -iXL-2)aVb + O (J^J ' 



(5.21a) 
(5.21b) 
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in which we pose 

atiy) = Xt 1 + {-fXt 1 . (5.22) 

Here we define X 1 = ^ = ±(1 + A) and X 2 = =a = 1(1 — A) with A = = ±J1 - Au, 

so that X l + X 2 = 1 and XiX 2 = z/. 21 

In addition we shall need the mass monopole / agreeing with its canonical counterpart 
M which parametrizes the various tail terms in Section IY Al Since the tails arise at 1.5PN 
order we need M only at the 1.5PN relative order. It is given by 

I = M = m(l-^)+o(^\. (5.23) 

We require also the current dipole moment or angular momentum Jj (agreeing with its 
canonical counterpart Si) since it appears in some non-linear terms, for instance in f)5.4p . It 
is needed only at Newtonian order, 

Ji = Si = vme iab x a v b + O \ K\ ■ ( 5 - 24 ) 

Finally, we have to provide the few gauge moments that enter the relations between 
canonical and source moments found in (I5.9P and (15.101) . They are readily computed from 
;he general expressions of all the gauge moments {Wl, X l , Y l , Z l } given in (5.15)-(5.20) of 



43j . The calculation is quite simple because these moments, namely the monopolar moment 



W and the two dipole moments Wi and Yj, are Newtonian. For circular orbits we find 

W = (J^J , (5.25a) 
Wi = ^z/mArV + ( — ) , (5.25b) 
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Note that the coefficient S(.(v) is equal to the product fhfk(v) in the notation of [19]. The equivalence of 



the two expressions follows from the Waring formulas [59J for X" + and X" — X%. We find (where 
(^j is the usual binomial coefficient) 



U 2fc-l (2k-l-p 
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or V cr 



We are done with all the source multipole moments needed to control the 3PN accurate 
FWF generated by compact binary sources in quasi-circular orbits. 

VI. TIME DERIVATIVES OF THE SOURCE MULTIPOLE MOMENTS 

For the purpose of computing the time derivatives of the source moments we require the 
3PN accurate equations of motion of compact binary sources. Like in the computation of 
the moments we have to take into account both the conservative effects at 1PN, 2PN and 
3PN orders, and the effect of radiation reaction at 2.5PN order. 

We consider non-spinning objects so the motion takes place in a fixed plane, say the x-y 
plane. The relative position x = yi — y 2 , velocity v = rfx/c/t, and acceleration a = d\/dt 
are given by 

x = r n , (6.1a) 
v = rn + rwA, (6.1b) 
a = (f — r uj 2 ) n + (r uj + 2r uj) A . (6.1c) 

For a while the time derivative will be denoted using an over dot. Here A = z x n is 
perpendicular to the unit vector z along the z-direction orthogonal to the orbital plane, and 
to the binary's separation direction n. The orbital frequency uj is related in the usual way 
to the orbital phase <fi by uj = <fi. 

Through 3PN order, it is possible to model the motion of the binary as a quasi- circular 
orbit decaying by the effect of radiation reaction at the 2.5PN order. This effect is computed 
by balancing the change in the orbital energy with the total energy flux radiated by the 
gravitational waves. At 2.5PN order this yields (see e.g. 



64 Gm r /9 _ / 1 i 
96 Gm 



uj = — ~^v~i 
5 r 6 



5/2 + O (i) , (6.2b) 
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where 7 is given by (15.131) . By substituting those expressions into (16. II) . 22 we obtain the 
expressions for the inspiral velocity and acceleration, 



a = —a; x 



5 
32 
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A central result of PN calculations of the equations of motion is the expression of the 



orbital frequency u in terms of the binary's separation r up to 3PN order. This resu 



been obtained in harmonic coordinates in 
in ADM coordinates in 
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58j and independently in 
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t has 



62], and 



65J. In the present work r is given in harmonic coordinates 



and the expression of the 3PN orbital frequency is 
r, Gm 
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Note that the logarithm at 3PN order involves the same constant r' as in the source 
quadrupole moment (I5.12p - fl5.14l) . This logarithm comes from a Hadamard self- field regu- 
larization scheme and its appearance is specific to harmonic coordinates. 

As often convenient we shall use in place of the parameter 7 given by (15.131) an alternative 
parameter x directly linked to the orbital frequency ( 16. 4p . namely 



x 



Gmw\ 2 ' 3 
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The interest in this parameter stems from its invariant meaning in a large class of coordinate 
systems including the harmonic and ADM coordinate systems. At 3PN order it is given in 
terms of x by 
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22 We notice that r = C(c 10 ) is of the order of the square of radiation-reaction effects and is therefore zero 
with this approximation. 
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Combining (I6.4p with (16.61) we find that the velocity squared v 2 = r 2 uo 2 + r 2 = r 2 uj 2 + 0(c 10 ) 
is related to x by 

= xi\ + x (-2 + ^) +x 2 (i + y^ + y) ( 6 - 7 ) 
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During the computation of the time derivatives of the source moments, each time an 
acceleration is produced the result is consistently order reduced, i.e. the acceleration is re- 
placed with ( I6.3bl) at the right PN order. Such an order reduction will generate in particular 
some 2.5PN radiation-reaction terms which are to be taken into account in the 3PN wave- 
form. This occurs when computing the time derivatives of the moments 7^, 1^ and that 
appear in the FWF at Newtonian and 0.5PN orders. On the other hand, when computing 
the polarization states following (12.21) we shall meet some scalar products of the polarization 
vectors P and Q with the relative velocity v. If those scalar products occur at Newtonian 
and 0.5PN orders {i.e. in multipolar pieces corresponding to the moments iy, Jyjj. and Jy) 
we shall have to take into account the 2.5PN radiation-reaction term coming from the ex- 
pression of v given by (I6.3a[) . 23 However it was shown in [3] that the radiation-reaction 
terms in the FWF at the 2.5PN order can be absorbed into a modification of the orbital 
phase, where they appear to constitute in fact a very small phase modulation, comparable 
with unknown contributions in the phase being at least of order 5PN - - negligible here 
since the phase is known only to 3.5PN order. In the present paper, we have chosen 24 to 
include all the radiation-reaction terms coming from both (I6.3ap and f!6.3bj) . and to present 



them as 2.5PN and 3PN amplitude corrections in our final results which will be presented 
in dSISD-dSUDD and (EHJ) below. 

Let us next check that the Hadamard self-field regularization constant r'o appearing both 
in the 3PN orbital frequency (16.41) and in the 3PN quadrupole moment (I5.14p . 25 is actually 



23 Not considering the radiation-reaction contribution in v given by (|6.3aj) has been the source of an error 
in [13] which has been pointed out and corrected in 0. 

As usual there are many different ways of presenting PN results at a given order of approximation, and 
choosing one or another is often a matter of convenience. 
25 The other moments are given at 2.5PN order at most; they do not depend on r'o since the appearance of 
regularization constants is a feature of the 3PN approximation. 
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a gauge constant. To this end we simply verify that r'o will be eliminated when expressing 
the FWF in terms of the gauge invariant parameter (16.51) . From (I5.14p we see that the 
dependence on r' of the 3PN quadrupole moment is 



Iij = vm 



44 n , f r 
l-- 7 3 ^ln — 
3 V. r o 



(6.8) 



We indicate by dots all the terms that are independent of r'o (for convenience we also 
show the Newtonian term). Now the FWF depends on the second time derivative of the 
quadrupole moment. For circular orbits this reads [coming back to the superscript notation 
(n) for time derivatives! 



l\f = 2um 



y7 3 v In (-^J (v (ij) + x (i a^ +--- + (J^J . 



(6.9) 



Replacing Vi and Oj by their values (16. 3p we get with the required approximation (still being 
interested only in the fate of the constant r'o) 

' 1 



I^f 1 = 2umv z 



44 o , / r 
1 _ y vkil — 

3 \r' 



x (ij) ~ n (H) H VO 



(6.10) 



The squared velocity v 2 = r 2 uj 2 + 0(c~ 10 ) appears in factor. It is now clear that replacing v 2 
by its expression in terms of the parameter x following (16. 7p . we produce another logarithmic 
term containing r'o, namely 



2 2 
V = C X 



44 , / r 

1 H x u\n\ — 

3 Vr'o 



+ ■■■ + - 



(6.11) 



which will cancel out the dependence of the quadrupole moment on r'o at 3PN order (using 
the fact that 7 can be replaced by a; in a small 3PN term). Thus, finally, 



if = 2u mc 2 x (A< y> - n m ) + ' ' ' + (J?) 



(6.12) 



is independent on r'o, which means that this constant cannot affect any physical result at 
the 3PN order. 



VII. COMPUTATION OF THE TAIL AND MEMORY INTEGRALS 



The results of Sections IVTTyTI yield the complete control of the instantaneous part of the 
FWF. We now tackle the computation of the hereditary part, which is composed of tails (and 
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tails-of-tails and squared-tails) and non-linear memory terms. The hereditary integrals have 
been explicitly provided in Section IV Al as contributions to the various radiative moments U L 
and Vl given by (I5.4l) - (l5.7p . Our computation will basically be a straightforward extension 
of the computation performed at 2.5PN order in Section 4 of [17 j. Since we employ exactly 
the same techniques, we skip most of the details and rely on 17] for justification of the 
method and proofs. 

We first consider the non-linear memory terms. Up to 3PN order we have the 2.5PN 
memory integrals in the radiative mass quadrupole moment Uij given by (15. 4p and the ra- 
diative mass hexadecapole moment Ujju given by (I5.6ap — these are the memory terms 



contributing to the FWF at 2.5PN order 



17| — and, in addition, we have the memory inte- 



gral in the mass octupole moment given by ( 15.5al) and the one in Uijkim given by ( I5.7af) 



- these contribute specifically at 3PN order. 26 Like in [17j we obtain the corresponding 
integrands (i.e. the terms under the integral sign) and compute directly their contributions 
to the two wave polarizations h + and h x . Indeed it is convenient to perform the relevant 
contractions of the integrands with the polarization vectors P and Q (see Section IVIIII for 
the conventions we adopt) so as to only deal with scalar quantities. 

We find that the memory integrals in h + and h x are composed of two types of terms. 
First there is a term, only present in the plus polarization h + , which does not depend on 
the orbital phase and can thus be viewed as a zero-frequency (DC) term. Actually, because 
of the steady inspiral, this term is a steadily varying function of time, with an amplitude 
increasing like some power law of the time remaining till the coalescence. Strictly speaking, 
this term is to be regarded as the memory contribution because it does depend on the 
behaviour of the system in the remote past, and therefore must be computed using some 
model for the evolution of the binary system in the past. In the present paper we find that 
the only zero- frequency term up to 3PN order is the one which appeared already at 2.5PN 
order and was evaluated in [17j — interestingly there are no other terms of this type at the 
3PN order. Because of the cumulative effect of integration over the whole past we know 
that this term, thou gh o riginating from 2.5PN order, finally contributes in the FWF at the 



Newtonian level 



5j, 



53 



541 ] . In practice the computation of this DC term reduces (in the 



Recall that the non-linear memory terms occur only in the mass-type radiative multipolc moments Ul- 
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circular orbit case) to the evaluation of the single elementary integral 

Here r(r) denotes the binary's separation at any time t <Tr (where Tr = T — R/c is the 
current time) . The coefficient in front of (17. ip is chosen for convenience to make the integral 
dimensionless. The integral (17.11) is easily computed using a simplified model of binary 
evolution in the past in which the orbit is assumed to remain circular apart from the gradual 
inspiral at any time. In this model the binary separation evolves like r(r) oc (T c — t) 1 / 4 where 



T r denotes the instant of coalescence (see 



171 ] for more details). In the remote past we thus 



have r(r) ~ (— r) 1//4 so the integral ( 17. II) converges when p > 4 (actually we shall only need 



the case p = 5 like 



inQ). 



The result reads 

where x(Tr) denotes the current value (i. e. at the current retarded time Tr) of the parameter 
x defined by (I6.5p . Witness the memory effect: the end result (17. 2p is of order x p ~ 4 = 
0(c~ 2p+8 ) which is a factor c 5 larger than the original formal PN order 0(c~ 2p+3 ) as shown 
in (17. ip . Hence, although the memory term is formally of order 2.5PN, its actual contribution 
to the waveform is comparable to a Newtonian term. As mentioned above we do not find 
memory (zero-frequency) contributions originating from the next 3PN order, and therefore 
finally no DC term at 0.5PN order. 

Second there are other terms, present in both polarizations, which depend on the or- 
bital phase, and oscillate like some harmonics of the orbital phase (say ncf)). Such phase- 
dependent, oscillating terms do not exhibit the memory effect, essentially because the oscil- 
lations, due to the sequence of orbital cycles in the entire life of the binary system, more or 
less compensate each other. As a result these terms, in contrast with (I7.1I) -( J7T21) . keep on 
their formal PN order. We recover the 2.5PN terms investigated in [17J and in addition we 
obtain several other terms at 3PN order. The latter are computed by a slight generalization 
of the method followed in jl^]: instead of (4.18) in [l^] we need to consider the integral 

XT K )= ( ^l_Jr — , (7.3) 

where 0(r) is the orbital phase at any time, where n and p range over integer or half-integer 
values {e.g. n = 1, 3, 5 and p = 11/2 at 3PN order), and where the coefficient is chosen to 
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make the integral dimensionless. Following the steps (4.18)— (4.23) in [l^] we compute this 
integral using our model of binary's past evolution, and in the adiabatic limit, which means 
that the current value of the adiabatic parameter £ associated with the binary inspiral is 
considered to be small and of PN order £(Tr) = 0(c~ 5 ). We then find 

Jn<f>(T R ) 



J{T R )=a*-i(T R )- 



in 



O 



(7.4) 



This result (valid only if n ^ 0) permits to handle all the phase-dependent oscillating terms 
coming from the memory integrals. 

We next turn to the computation of the tails and tails-of-tails present in the radiative 
moments (15.4p - fl5.7p . Again we closely follow the previous investigation [l^] on which we 
refer for more details. The computation of tails reduces to the evaluation of an elementary 
integral involving a logarithmic kernel, 

(Cm)'-' f T « { T r -t \ 
K0W = J dr In (^-^ j , (7.5) 



in which the logarithm has been scaled with the constant time T c — Tr, instead of the 
previous normalization by 2r , where T c is the instant of coalescence in the model of \^\ - 
Such scaling can always be done at the price of adding another term proportional to some 
integral of the type J(Tr) computed previously. Following the derivation of this integral in 
we find that, at dominant order in the adiabatic approximation, 



p in<t>(T R ) 

K(T R ) = x^(T R ) 



i n 



"-In fjLVc + /b "' 



(7.6) 



.2i \£(Tr) J V c 5 

Here C = 0.577 • • • is the Euler constant, and £(Tr) denotes the current value of the adiabatic 
parameter associated with the inspiral, which is defined by £(Tr) = \{T C — Tr)uj(Tr)}~ 1 in 
the model of \v\ . The adiabatic parameter is related to the PN parameter x by 

aT R )= 2 -^x^(T R ). (7.7) 
5 

The squared-tails are computed using the same integral (I7.5p - fl7.6p . Concerning the tails- 
of-tails we simply have to consider an integral involving a logarithm squared, 
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which is computed using the same technique with the result 



rP~9 



oincf>(T R ) 



in 



17T 

~2 



+ o 



In c 



(7.9) 



We are done with the computation of all tails and tails-of-tails in the 3PN waveform. 

For completeness let us give also the two technical formulas which enable ones to arrive 
at the results (17. 6p and (17. 9p . Posing y = (Tr — r)/{T c — Tr) and A = n/£, and working at 
the leading order in the adiabatic limit £ — > or equivalently when A — > +oo, the formulas 
express that, for any positive or negative A (see e.g. [66j p. 573 and 574), 



dy \iaye 



-i\y 



dy In 2 y e 



-iXy 



1 

A 

% 
A 



7T 



-sign(A) +z(ln|A| + C) 



° 1 i 



7T 

IF 



7T 



-sign(A) + i(ln |A| + C) 



(7.10a) 
(7.10b) 



Notice that we are only interested in the recent past contribution to the integrals (17.101) . 
corresponding to the interval < y < 1 equivalent to the time interval 2Tr — T c < t < Tr. 
The reason is that the remote past contribution, given by 1 < y < +oo or equivalently 
— oo < r < 2Tr — T c , is small in the adiabatic limit. This is a characteristic feature of tails: 
they die out very rapidly, therefore they depend essentially on the recent past evolution of 



the matter source 
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67|] . In the case at hand this technically means that the remote-past 



contributions to the integrals are of order 

dy In y e 



dy In ye 



-iXy 



-i\y 



0i V 



(7.11a) 
(7.11b) 



as can easily be verified by using integration by parts. 

VIII. 3PN POLARIZATION WAVEFORMS FOR DATA ANALYSIS 

We specify our conventions for the orbital phase and polarization vectors defining the 
polarization waveforms (12.21) in the case of quasi-circular binary systems of non-spinning 
compact objects. If the orbital plane is chosen to be the x-y plane (like in Section [VTl) . with 
the orbital phase <ft measuring the direction of the unit vector n = x/r along the relative 
separation vector, then 

n = x cos (ft + y sin , (8.1) 
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where x and y are the unit directions along x and y. Following [if], we choose the 
polarization vector P to lie along the x-axis and the observer to be in the y-z plane with 

N = s i y + c i z, (8.2) 

where we pose q = cosz and Sj = sini, with i being the orbit's inclination angle (0 < i < 7r). 
With this choice P lies along the intersection of the orbital plane with the plane of the sky in 
the direction of the ascending node A/", i.e. that point at which the bodies cross the plane of 
the sky moving toward the observer. The orbital phase <fi is the angle between the ascending 
node M and the direction of body one (say). The rotating orthonormal triad (n, A, z) 
describing the motion of the binary [see (16.11) ] is then related to the fixed polarization triad 
(N,P,Q)by 

n = P cos (f> + (q Q + Si N) sin <p , (8.3a) 
A= -P sin0+ (ciQ + SiN) cos0, (8.3b) 
z = - Si Q + c i N. (8.3c) 

As in previous works [if], we shall present the wave polarizations (I2.2p as expansion 
series in powers of the gauge-invariant PN parameter x defined by (16. 5ft . With a convenient 
overall factorization we write them as 




2Gmux 




4) 



c 2 R 

with the following PN expansion series 

ff+.x=E ajB/2fl Sx ) - (8-5) 

n=0 

The PN coefficients H { +^ ] will be given as functions of the orbital phase 0, and will also 
be polynomials in the symmetric mass ratio v and depend on the inclination angle i. In 
addition they will involve, at high PN order, the logarithm of x as we shall discuss below. 

Following [if], 3] it is convenient to perform a change of phase variable, from the actual 
orbital phase <fi satisfying <p = u, to some new variable denoted ip. Recall that the orbital 
phase evolves by gravitational radiation reaction and its expression as a function of time 
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is known from previous work 



13 



14, UM 

up to 3.5PN order. We then pose 



27 



2GMuo , / uo 
ip — <p 5 — m 

■ UJ 



(8.6) 



where M is the binary's total mass given by (15.231) . and where ujq denotes the constant 



ei2 



-c 



u 



4r 



•7) 



Here tq = tq/c is the normalization of logarithms in the tail integrals of the radiative 
moments (I5.4l) - (l5.7p ; r is the constant included in the definition of the finite part in (13.71) . 
Like r the constant u is arbitrary, because it is linked to the difference of origins of time 
in the far zone and in the near zone, see ( 13.161) . For instance we can choose uo = vr/ se i S mic 
where /seismic is the entry frequency of some ground-based interferometric detector. Using 
( I5.23P and the notation (16.51) the new phase variable reads 





In | 


f x \ 


1 X 




L 2 - 




\x ) 



where x = C^"^ ) 2 / 3 . 28 Our modified phase variable (18.61) - 



(8.8) 



will be valid up to 3PN 

n 



order but in fact it turns out to be the same as at the previous 2.5PN order [171 ] . 

The logarithmic term in ip corresponds to some spreading of the different frequency 
components of the wave along the line of sight from the source to the far-away detector, and 
expresses physically the tail effect as a small delay in the arrival time of gravitational waves. 
However, practically speaking, the main interest of this term is to minimize the occurence 
of logarithms in the FWF. Indeed we notice that the logarithmic term in (I8.6p . although 
of formal PN order 0(c~ 3 ), represents in fact a very small modulation of the orbital phase: 
compared with the dominant phase evolution whose order is that of the inverse of radiation 
reaction, i.e. <fi = 0(£ _1 ) = 0(c 5 ), this term is of order 0(c~ s ) namely 4PN in the phase 
evolution, which can be regarded as negligible to the present accuracy. Thus the logarithms 
associated with the phase modulation in (18.61) will be "eliminated" from the FWF at 3PN 
order. This does not mean that we should ignore them but that the formulation in terms of 
the small phase modulation (18. 6p is quite natural (for the data analysis it is probably better 



A similar phase variable is also introduced in black- hole perturbation theory [68j, |69|, |70( . 



We have lna; = § - |C - f ln2 



g In ( j§^- ) in agreement with the equation (68) of [19 1 
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to keep the logarithm as it stands in the definition of the phase variable ip). However all 
the logarithms will not be "removed" by this process, and we shall find that some "true" 
logarithms remain starting at the 3PN order. Such logarithms cannot be absorbed into some 
small modulation of the orbital phase, so 3PN will remain as the true order of magnitude 
of these logarithms in the FWF. 

With those conventions and notation we find for the plus polarization 29 
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29 We also requote the previous 2.5PN results [lTj taking into account the published Erratum [17| and the 
correcting term associated with radiation reaction and pointed out in [18( . 
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Notice the obvious fact that the polarization waveforms remain invariant when we rotate by 
7T the separation direction between the particles and simultaneously exchange the labels of 
the two particles, i.e. when we apply the transformation (ip, A) — ► (-0 + tt, —A). Moreover, 
due to the parity invariance, H + is unchanged after the replacement i — > tt — i, while H x 
being the projection of hj^ on a tensorial product of two vectors of inverse parity types, is 
changed into its opposite. 

We have performed two important tests on these expressions. First of all we have verified 
that the perturbative limit v — > of the polarization waveforms fl8.9p - fl8.10p is in full 
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agreement up to 3PN order with the result of black-hole perturbation theory as reported 



in the Appendix B of 69(. 30 Our second test is the verification that the wave polarizations 
f l8.9p - fl8.10p give back the correct energy flux at 3PN order. The asymptotic flux is given in 
terms of the polarizations by 



GW 



R 2 c 3 
hm — — 



dQ 

47T 



(h+) 2 + {h x y 



[8.11) 



where dfl is the solid angle element associated with the direction of propagation N. We 
have dQ = sinQdOdQ where (6,$) are the angles defining N, following the notation of 
Section [III To obtain the polarizations corresponding to this general convention for N we 
have to make some simple replacements in fl8.9l) - fl8.10l) for i and ip. As is clear from the 
geometry of the problem we must replace (i, ip) — > (©, xj) + tt/2 — <3>). The time derivative of 
the polarizations is computed in the adiabatic approximation, using <fi = uo and Co given by 
( 16.2bp . Of course one must take into account the difference between and the variable 
used in fl8.9l) -( l8TT0l) . Finally, the angular integration in (18. lip is readily performed and the 



result is in perfect agreement with the 3PN energy flux given by (12.9) of [l3j]. 31 

As already mentioned there are some "true" logarithms which remain in the FWF at the 
3PN order — i.e. after it has been expressed with the help of the PN parameter x and the 
phase variable ip. Inspection of (I8.9I) - (I8.10I) shows that these logarithms have the effect of 
correcting the Newtonian polarizations in the following way: 

1+ cf) cos2^j \ / 42g 



f H +\ 






H 




|05 xMn(16x)j+--- + (^-) , (8.12) 
-2Q sin 2ip x ' 

where the dots represent the terms independent of logarithms. In our previous computation 
of the 3PN flux using (18.111) we have already checked that these logarithms are consistent 
with similar logarithms occuring at 3PN order in the flux. Indeed we easily see that they 
correspond in the 3PN flux to the terms 



r 
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In a misprint was spotted in the Appendix B of the sign of the harmonic coefficient should 
be changed, so that one should read C 7 ^ 3 = + i 02 ?o 9 240 c °s(6 | )(167 + ...^inj^)^ 5 cos(3V>) - •••). 



The ambiguity parameters therein are now known to be A 



15] 



1987 
3080 
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65J and 9 = £ + 2k + ( 



11831 
9240 
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already known from (12.9) in 13J. Technically the logarithm in (18.121) or (I8.13P is due to the 



tails-of-tails at 3PN order. Notice that this logarithm survives in the test-mass limit v — > 



and is therefore also seen to appear in linear black hole perturbations 



68, 69, 



IX. 3PN SPHERICAL HARMONIC MODES FOR NUMERICAL RELATIVITY 

The spin-weighted spherical harmonic modes of the polarization waveforms at 3PN order 
can now be obtained from using the angular integration (I2.5p . An alternative route would be 
to use the relations (I2.6p - (I2.7I) giving the modes directly in terms of separate contributions 
of the radiative moments Ul and Vl- In the present paper the two routes are equivalent 
because all the radiative moments are "uniformly" given with the approximation that is 
necessary and sufficient to control the 3PN waveform. 

In this respect one should be careful about what we mean by controlling the modes up 
to 3PN order. We mean — having in mind the standard PN practice — that the accuracy 
of the modes is exactly the one which is needed to obtain the 3PN waveform. Thus the 
dominant mode h 22 will have full 3PN accuracy, but higher-order modes, which start at 
some higher PN order, will have a lower relative PN accuracy. For instance we shall see that 
the mode h M starts at 1PN order thus it will be given only with 2PN relative accuracy. 

The angular integration in (I2.5P is over the angles (G, $). Like in our previous computa- 
tion of the flux (18.111) . it should be performed after substituting (i, ip) — > (9, ip + tt/2 — $) 
in the wave polarizations. Denoting h = h + — ih x the integral we consider is thus 

h £m = j ' dn h(e, i) + tt/2 - $) Fl™(9, $) . (9.i) 

Changing $ into if) + tt/2 — if)' and into i' = arccos c' { , and using the known dependence of 
the spherical harmonics on the azimuthal angle $ [see (12. 4p ], we obtain 

h tm = (_^m e -im^ J ^ J ^ ^ ^ Y*™(i' , if)') , (9.2) 

exhibiting the azimuthal factor e~ m ^ appropriate for each mode. Let us factorize out in 
all the modes an overall coefficient including e~ im ^, and such that the dominant mode with 
(£,m) = (2,2) starts with one (by pure convention) at the Newtonian order. Remembering 
also our previous factorization in (18.41) we pose 

h em = 2G ™ c " X H em , (9.3a) 
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and list all the results in terms of H lm ^ 2 
107 55i/ 



H 



22 



1 + X 



+ x 



+ 



5/2 



42 

1077T 



+ 



42 



+ 27ra 3/2 + r 



2173 1069z/ 2047z/ 2 
+ 



1512 216 1512 
/ 27027409 856C 428ivr 



34ttz/\ , 

24w H I ~\~ oo i 

21 21 y V 646800 105 



+ 



105 



+ 



278185 41tt 2 



33264 



+ 



H 2L = -iA 



x 1/2 + ,r 



96 

3/2 



20261z/ 2 114635z/ 3 428 



2772 



+ 



99792 



105 



ln(16x) ) +0 



17 5i/ 

"28 + y 1 + r 1 :r + '' 



-- — 21n2 



5/2 



43 509z/ 79z/ 2 
+ 



17 

;1 56 +Z/ 



126 126 ' 16 
995 3 In 2 



H 20 = 
H 33 = 



14^6 



+ o 



84 
1 



+ 



x — 



17 In 2 



17"7T 37TZ/ 
+ 



28 



14 



14 



+ o 




3 .... . 
— i\ — A 

4 V 14 



x l/2 + x 3/2 ( _ 4 + 2z/) + ^2 + • + g ln (3/2 ^ 



= /123 1838z/ 887z/ : , 



127T + 



9ttu 



84 

/ | y -24 In (3/2) +u 



48103 
" 1215 



+ 91n(3/2) 



H 



32 



1 /5 



+ X 



3 V 7 

3 



x(l — 3v) + 



193 U5u 73z/ 2 



90 



+ 



1* 



+ o 



x 5/2 ( 2tt - 6ttz/ + 



1451 17387z/ 5557z/ 2 5341z/ 3 



3960 3960 



220 



1320 



+ o{-, 



H 31 = 



iA 



12V14 



a; 1 / 2 + x 3 / 2 



* ^W(ir + i(4-21n2 



5/2 



607 136i/ 247z/ 2 



198 



99 



19* 



'56 161n2 

-''T5 + — + !/ 



+ X 6 

I 71n2 

15 + 3 



87T 77TZ/ 
_ 



6 



+ o 



H 30 = -li\llx 5 / 2 v + o(\ 



9 V 7 



x(l — 3i/) + a; 2 



593 1273z/ 175z/ 2 



110 



+ 



66 



22 



32 



— im 



The modes having m < are easily deduced using # ■ m = {—) H 



51 



+ x 5/2 ( 4vr - 12nu + 



+ x* 



42 



1193 



241n2 +81n2 



1068671 10881191/ 146879z/ 2 226097z/ E 



43 



200200 
9i A 



28600 



2340 



17160 



(9.4h) 



4V70 



^ (1 _ 2 , )+;c V M _!£ + 1267, 131,= 



11 132 



33 



+ X 37T — 67TZ/ + i l-f 



32 /16301 



V 810 



- 12 In (3/2) j + 6 In (3/2) ) 



63 



+ i 



x(l — 3i/J + x 
21 84z/ 



437 805z/ 19z/ 5 



+ 



110 66 



22 



+ x 



5/2 



[ 2tt - 



67TZ/ 



(9.4i) 



x 



1038039 606751z/ 400453z/ 2 25783z/ 3 
+ 



200200 28600 



25740 



17160 



+ o 



H 



41 



84V10 



x 3/2 (1 _ 2z/)+a; 5/2 



+ x d 7T - 27TZ/ + i - — - 2 In 2 + v 



32 



15 



101 337z/ 83z/ 2 
~ 33~ 

+ 4 In 2 



33 44 
1661 



30 



+ o 



H 55 = 



504\/2 
625i A " 



(9.4j) 

(9.4k) 
(9.41) 



96V66 



x 3/2 (1 _ 2z/)+x 5/2 

181 



+ x* 5tt - IOttz/ + i 



14 



+ v 



263 688z/ 256z/ 2 
^9~ + ~39 39~ 
105834 



3125 



-20 In (5/2) +10 In (5/2) 



+ l^ 

(9.4m) 



H bA = 



32 



9V165 



(l -5z/ + 5z/ 2 ) + 



4451 3619z/ 521z/ 2 339^ 



+ 



# 53 = 



910 130 13 
69 464z/ 88i/ 2s 



+ 



26 



(9.4n) 



./jI A U (1 _ 2 „ )+I W_^ + ^_^ 

32 V 110 [ \ 13 39 39 J 



+ x \ 3n — Qnis + i 



543 / 83702 



V 3645 



70 



+ v 



12 In (3/2) +6 In (3/2) 



+ o - 



(9.4o) 



H 52 = 



H 



51 



27V55 
iA 



x 



(l - 5z/ + 5z/ 2 ) +x 3 



3911 3079z/ 413z/ 2 231z/ s 



910 130 13 
179 352z/ 4z/ 2 * 



+ 



26 



288 v / 385 

+ X 3 I 7T - 27TZ/ + % 



3/2 /-, on 5/2 / 179 352z/ 4l/ \ 

x 3/2 (l - 1v) +x 5/2 1 

V ; V 39 39 39 J 



181 

To" 



- 2 In 2 + v 



626 



+ 4 In 2 



+ o 



(9.4p) 



(9.4q) 



52 



£50 

H m 
H 65 

£64 

£63 
£62 
£61 
£60 

H 77 

£76 
£75 
£74 
£73 
£72 

H 71 

£70 
£88 
£87 

£86 



= O 



54 



5\/l43 
3125i x 5/2 A 



x 2 (1 - 5z/ + 5z/ 2 ) + x 3 



113 91z/ , P 39z/ 3 
+ — - 64z/ 2 + — — 



14 



128 A2" 



1 — 4i/ + 3i/ 2 



(9.4r) 

+ (i),(9.4 s) 
(9.4t) 



495 V 39 

81ia; 5/2 A 



x 2 (1 - 5i/ + 5z/ 2 ) + -— + 



93 71u tA 2 19z/ 3 
14 2 2 



616V65 

9 



1 — 4i/ + 3i/ 2 



+ 



297V65 

ix 5 / 2 A 



/ 81 59z/ 9 7z/ 3 

x 2 (1 - 5i/ + 5z/ 2 ) + x 3 ( -— + — - 32z/ 2 + — 



= O 



8316V26 
1 



1 — 4i/ + 3i/ 2 



+ 



16807ia; 5 / 2 A / 7 



1440 



858 



1 - 4i/ + 3i/ 2 



+ 



35 V 143 V \c 7 / 



15625i x 5/2 A 



26208V66 



1 - 4i/ + 3i/ 2 



128x 3 AT , P ^ / 1 

243i:r 5 / 2 A fs 



160160 



x 3 (1 - 7u + 14z/ 2 - Iv 



1 - 4i/ + 3z/ 2 

3^ 



3003^ 



+ K 



864864\/2 
16384 



1 - 4z/ + 3z/ 2 



63 
1 

7 



243 



-x 



35 V 17017 



(1 - Iv + 14z/ 2 - 7z/ 3 ) + ^ , 



(9.4u) 
(9.4v) 

(9.4w) 
(9.4x) 
(9.4y) 
(9.4z) 
(9.4aa) 
(9.4bb) 
(9.4cc) 
(9.4dd) 
(9.4ee) 
(9.4ff) 
(9.4gg) 
(9.4hh) 
(9.4ii) 
(9.4jj) 



53 



£85 

£83 
£82 
£81 
£80 



o 



o 



128 
4095 
1 



x: 



9009V85 



(l-7z/+14z/ 2 -7z/ 3 ) +0 , 



(9.4kk) 
(9.411) 
(9.4mm) 
(9.4nn) 
(9.4oo) 
(9.4pp) 



while all the higher-order modes fall into the PN remainder and are negligible. However, we 
shall give here for the reader's convenience their leading order expressions for non zero m 
(see the derivation in IS]). For £ + m even we find: 
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where we recall that the function s^(y) is defined in (15.221) . For £ + m odd we have: 

f5(£ + 2){2£+l){£ + m)\{£-m)\\ 1/2 



H 



2^-i ^+ m ~ 1 ' 



\(2£+ 1)! 



£(£-!)(£ + !) 



X 8 t+ i 



'y)i[im) l x {i - x)t2 + 



1 



(9.6) 



When m — 0, -ff^ m may not vanish due to DC contributions of the memory integrals. We 
already know that such an effect arises at Newtonian order [see (I8.9al) ]. hence the non zero 
values of H 20 and if 40 . 

We find that the result for H 22 at 3PN order given by (I9.4al) is in complete agreement 
with the result of Kidder 19] . The only difference is our use of the particular phase variable 
( 18. 8p which permits to remove most of the logarithmic terms, showing that they are actually 
negligible modulations of the orbital phase. For the other harmonics we find agreement with 
the results of [19| up to 2.5PN order, but the results have here been completed by all the 



3PN contributions. 
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APPENDIX A: LIST OF SYMBOLS 

a relative acceleration of binary masses in harmonic coordinates 

a im STF tensor connecting the usual spherical harmonics basis Y £m 

to the STF tensors basis 
Ci cos i 

A mass difference ratio; A = (mi — m 2 )/m 

JF GW total gravitational wave energy flux 

FWF full gravitational waveform at 3PN order 
7 PN parameter; 7 = 

hj? transverse traceless (TT) projection of metric deviation; Eq. (12. ip 

h +jX "plus" and "cross" polarisation states of the FWF; Eqs. (12.21) 

iJ +>x same as h +iX modulo an overall factor; Eqs. (18.51) . (18.91) and (18.101) 

h £m spin-weighted spherical harmonic modes of the FWF; Eq. (12.51) 

H £m , H £m same as h tm modulo overall factors; Eqs. ( j9.3f) 

Il mass-type source multipole moment STF with I multipolar spatial indices; 

is given for i = 2, 3, 4, 5, 6, 7, 8 by Eqs. (15TT2]) . (05aj] . 
(PngaD , d5TTa| . (JEHU), flSHID, respectively 
i inclination angle of the binary orbit 

Jl current-type source multipole moment STF with i multipolar spatial indices; 

is given for t = 2, 3, 4, 5, 6, 7 by Eqs. fj5T5bjl . (j516bl) . (Q7bjl . 

fl5~T8bj) . fl5~T9bl) . fl5T20bD respectively 
I multipolar order 

A unit vector in the orbital plane; A = z x n 
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mi , 1712 individual masses of binary components 

m total mass of the binary; m = mi + m 2 

M source mass- type monopole moment; Eq. (15.231) 

Ml canonical mass-type STF moment with I multipolar spatial indices, related to 

source moment II by Eqs. ( 14.25al) and ( 15. 91) . ( 15. 10a[) 
n binary's separation direction, from 7712 to m\ 

v symmetric mass-ratio; v = mi ^ 2 

N direction of propagation of gravitational wave 

P unit vector along the direction of the ascending node M 

Vj^. x TT projection operator; Eq. (12. ip 

$ azimuthal angle of N in spherical polar coordinates 

4>(t) orbital phase of the binary's relative orbit, the angle between n and P, 

increasing in the direction of A 
ip(t) effective orbital phase of the binary's relative orbit, as modified by 

4PN log term; Eq. (KTSjl 
Q unit polarization vector in the plane of the sky; Q = N x P 

RWF restricted post-Newtonian gravitational waveform 

r relative separation of binary masses in harmonic coordinates 

R distance to the source in radiative coordinates 

STF symmetric-trace-free projection 

Si sin i 

Sl canonical current-type STF moment with I multipolar spatial indices, related to 

source moment Jl by Eqs. (I4.25bj) and (I5.10bj) 
O polar angle of N in spherical polar coordinates 

TT transverse traceless projection 

Ul mass-type radiative multipole moment STF with i multipolar spatial indices; 

is given for £ = 2, 3, 4, 5 by Eqs. (15.4)) . (I5.5a|) . (I5.6al) . (I5.7aj) respectively 
U £m radiative mass moment (non-STF form) corresponding to h tm \ Eq. fl2.7|) 

v relative velocity of binary masses in harmonic coordinates 

Vl mass-type radiative multipole moment STF with £ multipolar spatial indices; 
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is given for i = 2, 3, 4 by Eqs. (jObl) . (l5Jb]) . fl5T7b|) respectively 

V^ m radiative current moment (non-STF form) corresponding to h £m ; Eq. (12.71) 

u angular velocity of the relative orbit in harmonic coordinates; Eq. ( 16. 4p 

Wl gauge moment, entering the relation between canonical and source moments; 

Eqs. (EH}, (EH 

X L gauge moment 

x gauge invariant PN expansion parameter; Eq. (16.51) 

x^ harmonic coordinate system in the near- zone 

X M radiative-type coordinate system in the far-zone 

yi™(0,$) gpin.-weighted spherical harmonics of weight —2; Eq. (12. 4p 

y £m (0,$) standard spherical harmonics 

Y L gauge moment 

Z L gauge moment 

z unit vector normal to the binary orbital plane 
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